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Abstract 

Following the construction of the general effective theory for dark matter direct detection in 
PQ , we perform an analysis of the experimental constraints on the full parameter space of elastically 
scattering dark matter. We review the prescription for calculating event rates in the general effective 
theory and discuss the sensitivity of various experiments to additional nuclear responses beyond 
the spin-independent (SI) and spin-dependent (SD) couplings: an angular-momentum-dependent 
(LD) and spin-and-angular-momentum-dependent (LSD) response, as well as a distinction between 
transverse and longitudinal spin-dependent responses. We consider the effect of interference between 
different operators and in particular look at directions in parameter space where such cancellations 
lead to holes in the sensitivity of individual experiments. We explore the complementarity of different 
experiments by looking at the improvement of bounds when experiments are combined. Finally, our 
scan through parameter space shows that within the assumptions on models and on the experiments' 
sensitivity that we make, no elastically scattering dark matter explanation of DAMA is consistent 
with all other experiments at 90%, though we find points in parameter space that are ruled out only 
by about a factor of 2 in the cross-section. 



1 Introduction 

The idea that the dark matter (DM) component of our universe is a new particle, or perhaps 
a new of set of particles, remains very attractive. Current and future direct detection exper- 
iments offer the exciting possibility of obtaining the most complete picture of the properties 
of this particle. Traditionally, the DM particle has been viewed mostly through the WIMP 
framework, whereby the DM particle is intimately connected with theories of electroweak 
breaking. Thus, theories of DM have been written in terms of relativistic degrees of freedom 
relevant for TeV scale physics. On the other hand, the DM detection experiments measure 
the signal due to a non-relativistic collision of DM with nuclei in the body of the detector. 
When comparing experiments to each other, it is therefore more prudent to use a description 
that does not suffer from model-building biases coming from assumptions about physics at a 
much higher energy scales. Indeed, such assumptions have often lead to a view that there are 



only two possibilities for the interaction between nuclei and DM. Namely, that the coupling 
is either spin-independent (SI) or spin dependent (SD). 

In a previous paper, extending the work of [2], we argued that with minor assumptions, 
there are, in fact, several non-relativistic operators consistent with the rules of QM, Galilean- 
invariance, and CP symmetry, which describe possible nucleon-DM interactions. An im- 
portant difference with the standard paradigm is that we allowed for momentum dependent 
interactions as well (to second order in the momentum). Such interactions are well- motivated, 
and fully relativistic models for these have been studied earlier as a way to reconcile the vari- 
ous experiments with each other. For example, the DM particle could be composite (after all, 
most visible matter takes the form of atoms), and its interaction with nuclei could proceed 
through dipole or charge-radius type operators sensitive to its size. The advantage of this 
framework is that it is essentially model-independent and makes use only of the ingredients 
that DM experiments are probing. The operators we focus on also apply for DM particles of 
any spin. Any high energy relativistic DM model can be translated into a particular linear 
combination of the EFT operators. On the other hand, the non-relativistic operators form 
a general basis for the interactions most directly probed by the experiments themselves. In 
many ways, these operators are analogous to the precision electroweak observables, which 
were useful in parameterizing the results of precision experiments measuring the properties 
of the W and Z in a model independent way. Thus, instead of thinking of target nuclei in 
terms of SI and SD features alone, a more model independent analysis should be one which 
considers the nuclear responses coming from the non-relativistic interactions. 

Interestingly, we found that some of the new momentum-dependent interactions give rise 
to novel nuclear responses not considered in the DM literature. In particular, there are five 
different response functions which enter. They are: the standard SI response, two types of 
spin- dependent responses, SD1 and SD2 (a particular combination of which is the usual SD), 
a response which for certain elements is significantly dependent on the angular-momentum of 
unpaired nucleons within the nucleus (LD), and a novel response which depends on a product 
of spin and angular momentum (LSD). In addition, interference is possible between operators 
with similar quantum numbers. In certain regions of parameter space, this interference could 
significantly alter the strength with which DM couples to a particular target nucleus. 

In this paper we will first present the non-relativistic operators and the five nuclear re- 
sponses they elicit. Our goal will then be to explore the parameter space of these operators, 
given current experimental bounds. In particular, as various targets will couple differently to 
DM, our focus will be to determine which experiments need to be combined in order to opti- 
mize sensitivity to a given DM scenario. It is important to check whether there are potential 
"blind spots" with the targets currently being used, In addition, it is useful to know which 
targets compliment each other the most, yielding better sensitivity in the future. 

Our paper is organized as follows. In Sec. [2] we review the non-relativistic effective theory, 
listing the operators and present our framework for analysis of the responses. In Sec. [3] 
we present the current experimental constraints on the parameter space, emphasizing where 
there are potential gaps in sensitivity. We also provide a qualitative description of the various 
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responses. Finally, in Sec. [4] we summarize our results and conclude. 

2 Review of Effective Theory 

2.1 Building Blocks of the Theory 

Interactions in the effective theory are constructed from a small number of "building blocks" 
that are directly connected to the non-relativistic process of a dark matter particle scattering 
off of a nucleon inside a nucleus. The first of these building blocks are the momenta of the 
particles in the collision. Since we are interested in two-to-two scattering, there are a priori 
four different momentum vectors, namely p, k for the incoming dark matter and nucleon 
momenta, respectively, and p', k' for their outgoing momenta. However, the constraints of 
momentum-conservation and independence of inertial frames allows us to reduce these from 
four to two independent momenta. Without loss of generality, one may take these to be 
q = p ' — p, which is the momentum transfer of the collision, and v = v x ^ n — VN,m, which is 
the incoming DM velocity in the nucleon rest frame. While these are both manifestly frame 
independent quantities, they are not the most convenient combinations to consider, because 
the most general interaction built out of them is not Hermitian. Since Hermitian conjugation 
interchanges in and out states, one effectively has q — > —q and v — >■ v y out — out = v + — , 
where fi^ is the dark- matter- nucleon reduced mass. Consequently, it is natural to work 
instead with the quantities 

iq and v 1 " = v H . (1) 

2/iTV 

In addition, the spins S x and Sn of the dark matter and the nucleons, respectively, are 
Hermitian operators that can appear in interactions. For spin-| fermions (i.e. for the dark 
matter in many models, and the nucleons in any case), these are just Pauli sigma matrices. 
The most general set of interactions for elastically scattering dark matter are the rotational 
invariants that can be constructed out of these four basic building blocks: 

S x ,S N ,iq, and v L . (2) 

2.2 List of Effective Operators 

In practice, we will not study all possible rotationally-invariance combinations of the building 
blocks in eq. (|2]), but instead will make the further restriction that they can arise from the 
exchange of a spin-0 or spin-1 field. Effectively, this means including at most two powers 
of the spins and/or velocity. This leaves us with the following T-even operators, classified 
according to symmetry into different groups such that operators in different groups will not 
interfere with each other: 
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1. P-even, S^- independent 

Ox = l, 2 = (^) 2 , 3 =iS N -(qxv x ), (3) 

2. P-even, ^-dependent 

0, = S x -S N , 5 = ^4-(gx^), 6 = (4-g)(^-g), (4) 

3. P-odd, 5 X - independent 

O r = S N -v ± , (5) 

4. P-odd, ^-dependent 

8 = £ x • tf 1 , 9 = 4-(^x 9 1 (6) 

In addition, there are T-odd operators [I], but since these do not introduce any new nuclear 
responses we will not analyze their constraints separately. 

2.3 Direct Detection Event Rates 

The calculation of direct detection rates in the general effective theory follows along essentially 
the same lines as in the case of standard spin-independent and spin-dependent interactions, 
which are reviewed in e.g. |3j. However, we will see that there are a few non-trivial gen- 
eralizations, and ultimately the final result will need to be parameterized differently from 
usual. For instance, it is common to parameterize direct detection rates in terms of the total 
wimp-proton cross-section, but more generally this quantity is momentum-dependent. The 
detector event rate Rd per unit time per unit energy per unit detector for a nucleus to recoil 
with energy Er is given by an average over the underlying dark matter velocity distribution 

dRo at / d& 



where n x is the local DM density, v x is its velocity in the lab frame, Nj< is the number of 
nuclei per detector mass, and a? is the scattering cross-section for dark matter and atomic 
nuclei. The kinematics of scattering are easiest to follow in the center-of-mass frame, where 
the dark matter momentum p and target momentum fc^ are equal and opposite, and given by 

p = —k,T = HtV Xi (incoming) (8) 

p' = -k' T = q + fi T v x , (outgoing) (9) 

where q = p'—pis the momentum transfer, and fiT is the dark matter- nucleus reduced mass 

,T = (10) 

ttit + rn x 
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In this frame, energy conservation implies \p\ = \p'\, so 
a 2 

y = p 2 -p-p' = p 2 {l - cos 9) = ii 2 T v 2 x (l - cos 9), (11) 

where cos 9 is given by the scattering angle in the center-of-mass frame. Then, the measured 
recoil energy is just the kinetic energy transferred to the nucleus: 

The minimum velocity required for a given momentum transfer occurs at exact backscattering 
(cos# = -1): 

^min = 7^—- ( 13 ) 

Thus, holding E R fixed and integrating the velocity over a halo distribution f(v) implies the 
following formula for the event rate: 

dR D Ar p x m T f o da T 

— — = Nt— — ^- / rf 3 w 1 f(v)- 1 - n . 14 

d^fl m x ^J Vmin rfcos^ 

The differential two-to-two cross-section is given as usual in terms of the matrix elements by 
averaging over initial spins and summing over final spins: 

d*T _J 1 V 1 \M\ 2 

d cos 9 2j x + 1 2j + 1 s ^ 32tt (m x + m T ) 2 ' 1 J 

where j and j x are the spin of the nucleus and dark matter, respectively. The calculation 
of these matrix elements depends on the ground state of the atomic nucleus, but the result 
can be parameterized in terms of form factors for the various effective theory operators. We 
parameterize the resulting amplitude by 

^^TEIMI 2 - ^/oFoM (16) 



spins 



The factor m^/m 2 N is explicitly factored out of the definition of Fq for convenience, due to the 
difference in relativistic normalization of states for nuclei (k T \k' T ) = (2ir) 3 (2m T )S( 3 \k T — k' T ) 
vs. nucleons (k\k') = (2ir) 3 (2m N )5( 3 \k — k'). More generally, in the presence of multiple 
effective interactions, there may be interference terms, each with their own form factor: 

^Vt^tEixi 2 . <j: CiCjF , j(vW) . (17) 

A spins JV ij 

Thus, putting everything together, the generalized event rate Rd per unit time per recoil 
energy per detector mass is 

dR D p x m T 



dE 



» = ^s^L*™- 1 ™^™™- (i8 » 
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2.4 A worked example: Dark magnetic moment 

As a practical example of the use of the effective theory and the form factors described in 
the previous section, we will now discuss in detail how to treat the case where dark matter 
couples to the Standard Model through a dark magnetic moment. That is, consider a dark 
sector that contains a massive gauge field A'^ that mixes kinetically with the photon 

C D eF'^, (19) 

and furthermore, the leading interaction of the dark matter \ with A'^ is through a magnetic 
moment interaction. The standard way of writing the magnetic moment interaction non- 
relativistically is as a coupling between the \ spin S x and the dark magnetic field B' = V x A: 

£D2m x /i DM ^-5 / . (20) 

This interaction by itself is not frame- in dependent, and must combine with other interactions 
to form a boost-invariant combination. When x is spin-^, such a combination arises from the 
following Lorentz-invariant operator: 

Ant = fXDMX^xKv ( 21 ) 

Integrating out the massive A' one obtains the interaction between \ an d Standard Model 
matter: 



CD C -^(xia^q aX m M , (22) 

where j£ M is the electromagnetic current. Restricted to protons and neutrons, it can be 
written 

J ™ = p{k } {-^r + -2^) p{k) +n{k) {-2^) n ^ (23) 

where g p = 5.59 and g n = —3.83 are the proton and neutron magnetic g-factors, respectively. 
Since coupling through kinetic mixing with the photon is a particularly compelling class of 
models, and this always leads to dark matter interactions through the EM current, we discuss 



these more generally and in more detail in App. [Bj Continuing with the interaction eq. (22), 
we can take its non-relativistic limit to obtain 



m 4 



proton 



m N q 2 l + 4m N m x iS x ■(qxv ± ) + 2g p m x (q 2 S x ■ S p -(q- S x )(q • S p ))) 

' i 

+ (2g n m x (q 2 S x ■ S n - (<f ■ S x )(q ■ S n ))) ] . (24) 

V / neutron. 

The first line are interactions with the proton, and the second are with the neutron. Com- 
paring with the definitions of the effective theory operators in Sec. 2.2, we can read off the 
coefficients for this model: 
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protons : = m N q 2 G M , cf ] = 2g p m x q 2 G M , c£ ] = 4m N m x G M , = -2g p m x G M , 

neutrons : C4 = 2g n m x q 2 GM, = —2g n m x GM, (25) 



where 



e/x DM e 

G M = - 2 ■ {2b) 



m A 



Summing over the coefficient, and using the formulae in App. [C], we find that 

£ aCjFvtf, q 2 ) = q'm\G M \m% ( -L + % _ ±) F ^\q 2 ) + AF^\q 2 ) (27) 
ij L V m x 5 /~r/ 

"2 (^<?(g 2 ) + 9 P F^)) + \ 9N9N^ N \q 2 ) 

N,N'=n,p 

3 Constraints and Gaps 

3.1 Qualitative discussion of Responses and Gaps 

We would like to first review the five responses of interest and provide a qualitative picture 
of the sensitivities of the various elements to each one. Due to interference, it will ultimately 
be useful to group the responses into sectors consisting of operators that can interfere with 
each other. However, let us first start by describing the novel responses individually. 



3.1.1 Standard spin-independent response (SI) 

The standard SI response, or M pn , is unchanged in our framework, and arises from 0\. A 
convenient measure for discussing the strength of various operators is just the interaction 
itself squared and averaged (summed) over initial (final) states, which we will label (£f nt )e>. 
For Ox, this is 

(A ? nt)oi ~ K N, (28) 
where is the coherence factor, defined as (A — Z),Z for N = n,p. 



3.1.2 Spin dependent responses (SD1 and SD2) 

Even in the case of spin-dependent interactions, there are in fact two different ways that the 
nucleus can couple to spin. The first one, SD1 (or S^ n ), is a projection of spin in the direction 
parallel to the momentum transfer, q, as in the case of operator Oq. The second, SD2 (or 
Ep n ), is a projection of the spin in a direction perpendicular to the momentum transfer, as 
for operator Og. Qualitatively, both these responses behave similarly to the standard SD 
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response, favoring elements with unpaired nucleons in the outer shell (p or n). Thus, C 6 and 
Og are roughly of the same relative size for different elements, except that they differ in their 
momentum-dependence: 

(C 2 Qt ) 0(i ~ q 4 S 2 N , (29) 
(AntK ~ <fSl- (30) 

However, quantitatively, SD1 and SD2 are different from each other and from the usual SD 
response (which arises from O4 = Sn • S x and is a linear combination of SD1 and SD2). In 
addition, as we will describe later, the SD2 response can interfere with responses sensitive to 
the angular momentum content of the nucleus (LD), producing important effects, which can 
reduce the sensitivity of targets to DM in various regions of parameter space. 



3.1.3 Angular-momentum dependent response (LD) 

Operators O5 and 0% contain the standard SI response M pn . However, for nuclei which 
contain an unpaired nucleon with angular momentum, there can be an important correction 
due to an angular-momentum-dependent response LD (or A p n ). This correction is of order 
one for elements with unpaired protons, such as 23 Na and 127 I, and for the isotopes with 
unpaired neutrons such as 73 Ge and 131 Xe. For C 5 , the approximate interaction strength is 

(Oft - ^( L n + K^). (31) 

2 

Note that for m x > m T , the ratio ^ is approximately l/A 2 , so that both terms are roughly 
comparable. Also, of qualitative importance, is the fact that the LD response can interfere 
with SD2, reducing the sensitivity to DM in certain regions of parameter space. 

3.1.4 Angular momentum and spin dependent response (LSD) 

The dominant response for operator O3 is one which is sensitive to a nuclear feature that has 
not been considered previously, namely the product of spin and angular momentum (L ■ S) p>n : 

4 ^2 
(A ? nt)o 3 ~ ~It(( l n ■ S N ) 2 + S 2 N —f). (32) 



Thus we will refer to it as the LSD response (or $p„). As above, ^ tends to suppress the 
second term, so the LSD response tends to dominate. To get a sense for this response, recall 
that when all 2(£ + 1) states of the spin- aligned (j — £ + \) subshell and all 21 states of the 
spin- ant i- aligned (j — £ — |) subshell are occupied, this dot product vanishes. However, the 
I = j + \ orbital and t — j — \ have different energies, so one will start to fill before the 
other. Taking n±(£) to be the occupation numbers of the £± \ orbitals, the expectation value 
(Sn-L n ) is proportional to (£ + 1) n + (£) — £ n_(£). Thus a mismatch between n± (£) of order 
£ produces ((L ■ S)) ~ ^highest f° r m °st elements. £ 2 hi g hest grows with the atomic number of 
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the nucleus, and so this response tends to favor heavier elements. Qualitatively, it is thus 
somewhat similar to the standard SI response, although important differences occur. For 
example, 19 F is far less responsive than 23 Na, even though their SI properties are similar. An 
additional important qualitative aspect of the LSD response is that it can interfere readily 
with the standard SI response, should the latter have non trivial momentum dependence (as 
can be natural in UV relativistic models which generate O3 in the first place). 



3.2 Constraints on different responses 

Because the nuclear responses to the full set of effective theory operators in eqs. (|3])-([6]) 
depend only only a smaller set of 5 independent responses SI, SD1, SD2, LD, and LSD (i.e. 
M, A, and $"), we will begin by showing the constraints on a representative effective 

theory operator for each of these. These constraints will have overlap with those considered 
in [2], for which operators involving SI, SD1 and SD2 were considered. Fig. [I] therefore shows 
constraints from individual experiments on Oi, O3, O5, Oq, and Og. By consulting eqs. (50a)- 
(501), one can read off that these five effective theory operators are sensitive to SI, {LSD 
and SD2}, {SI and LD}, SD1, and SD2, respectively. Note that no operator is sensitive to 
LSD alone or LD alone - these always appear with SD2 or SI, respectively. The following 
salient features emerge, in agreement with our above qualitative discussion of the different 
types of responses. First, the strongest constraints in most cases is from XENON100, mainly 
due to its significant exposure and small background. The only exceptions are C 6 and Og 
acting on protons, which are sensitive only to the relatively small proton spin in xenon-129 
and xenon-131. To separate out the effect of the large exposure vs. the intrinsic sensitivity 
of the isotopes to various operators, note the following experimental approximate exposures 
after including acceptances and efficiencies: 

effective 
exposure (kg d) 
CDMS I ~ 200 

XENON100 " ~ 2500 



COUPP 



DAMA 



25 



50/q 



In the case of DAMA, the sensitivity is determined not by total exposure but by the amplitude 
of the irreducible modulating component that they observe, which depends on the quenching 
factor q. Depending on whether the modulating component is due to scattering off of sodium 
or iodine, the appropriate quenching factor is gNa ~ 0.3 or qj ~ 0.085, respectively. Similarly, 
COUPP sees 20 nuclear recoil events [I], which are all treated as potential dark matter events, 
so we choose their effective exposure in the above table to be their total exposure divided 
by 20. To obtain an estimate for the constraints, one can estimate the number of predicted 
events as 



dN 
dE~ R 



5000keV 



.! / exposure \ /100GeV\ 
V kg • day J \ m x J 



(33) 
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For the high-energy analyses, the the effective exposures above have all been defined so that 
the constraint is roughly j^- < 1. 

In p], various nuclear models were used to calculate the relevant nuclear form factors; we 
can read off from them the size of the interaction terms above: 





Q2 
D n 


C2 
D p 


T 2 


T 2 


(S n ' L n ) 2 


(S p ■ L p ) 2 


F 


8 • 10" 5 


0.2 


0.04 


0.05 


0.6 


0.1 


Na 


0.0004 


0.06 


0.1 


0.8 


5.5 


3.3 


Ge 


0.02 


5 ■ 10" 6 


1.1 


0.003 


35 


100 


I 


0.004 


0.07 


0.4 


2. 


100 


500 


Xe 


0.02 


2 • 10~ 5 


0.4 


0.04 


500 


300 



All isotopes have been averaged over according to their natural abundance. For compari- 
son and to give a sense of the size of the uncertainties involved in the calculation at small 
momentum transfer, the following table provides a summary of the quantities S 2 , S 2 , L 2 , L 2 n 
in various treatments in the literature. Reference [5] present the results from two nuclear 
models, both presented below. The differences at fluorine and sodium are negligible, whereas 
the differences in the heavier elements - especially with protons for germanium and xenon 
and with neutrons for iodine - can be quite significant. 





C2 
°n 


C2 

o p 


T 2 


T 2 
L P 


F| 


8 -10" 5 


0.2 


0.04 


0.05 


Na 


0.0004 


0.06 


0.1 


0.08 


Ge[7! 


0.02 


9-10" 6 


1.0 


0.02 


I [5] i) 


0.006 


0.09 


0.6 


2. 


ii) 


0.004 


0.1 


0.4 


2. 


Xe [5J i) 


0.04 


0.0002 


0.5 


0.02 


ii) 


0.03 


7 -10" 5 


0.5 


0.05 



3.3 Interference sectors 

We will now turn to exploring the full parameter space of interactions, where multiple in- 
teractions can contribute to scattering. As a result of interference between operators, when 
considering current experimental constraints, it is most useful to quote these constraints in 
terms of a bound on the sum of coefficients squared (i.e. Yli c l) f° r each interfering sector. In 
other words, we will be looking for flat directions within each sector, trying to determine the 
size of the largest coefficients allowed in front of interfering operators of identical dimension. 
We will restrict our analysis to the case where the operator coefficients in the effective theory 
are constants, independent of the momentum-transfer. In principle, one could add additional 
q 2 dependence to the coefficients and parameterize the result as in [8]. But we do not treat 
this case here. 
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0[ Nucleon Scattering 




O5 Nucleon Scattering Of, Nucleon Scattering 

Log 10 (c 2 /GeV- 2 ) Log 10 (c 2 /GeV- 4 ) 




Figure 1: Constraints on representative operators of the five independent nuclear responses, 
for each individual experiment. 
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One class of cancellations, coming from the interference in isospin-dependent interactions, 
has been considered in attempts to explain the current picture of direct detections (see e.g. 
[9]). In a model-independent framework, isospin is simply added to the parameter space 
of operators as a separate proton and neutron interaction coefficient, c p and c n , for each 
operator. However, there is additional interference that occurs even between operators of 
the same isospin. In particular, the LD response interferes with the SD2 response (but both 
have the wrong parity to interfere with SD1). The interference occurs entirely due to the fact 
that nucleon-DM operators containing the velocity operator v are sensitive to the internal 
motion of nucleons within nuclei. Roughly, one can think of this as due to the fact that for 
the relative motion of a nucleon, N, around the nucleus, VN,rei ~ ^r~9 x ^n, where is the 
angular momentum of that nucleon. Additionally, the SI response interferes with the LSD 
response, because the expectation ((L • S)) is non-zero for most nuclei. 

In practice, the constraints of course depend not only on the intrinsic sensitivity of the 
nucleus to the interaction but also on experimental setups such as exposures and efficiencies; 
a large exposure can compensate for a small intrinsic interaction strength. Calculating the 
constraints on the full parameter space will address the question of whether current experi- 
ments are complementary in probing the full range of interactions. We will also discuss in this 
section the importance of having a diversity of nuclear targets. An important consideration 
in understanding constraints in the full parameter space is the dimensionality of interactions 
that an isotope can in principle be sensitive to. This dimensionality is limited by the spin 
j of the nucleus, since the maximum partial wave of a response that can contribute is 2j. 
Thus, if this number is less than the dimension of the parameter space, we can always find 
a "flat" direction along which the target is not sensitive due to interference. In general this 
direction varies with q and will be "smeared out" after we integrate over recoil energy. For 
light elements like F, however, since \q\ is small, the event rate is populated only at low recoil 
energy, and thus the flat direction remains. This explains the fact that in the interference 
sectors (Fig. [2]) there is no constraint coming from SIMPLE. In this case, another experiment 
that is sensitive in this direction will significantly improve the overall constraint. 

Similarly, we see the advantage of having many isotopes in the natural abundance of an 
experiment. This is because the isotopes will differ in their sensitivity to the various neutron 
interactions, with different sensitive directions. Generically, the flat direction of one isotope 
will be more sensitively probed by another, and hence the overall sensitivity for the element 
will be better than the individual isotopes. Thus, the experimental constraint will be stronger 
than the constraint coming from each individual isotopic fraction. Furthermore, isotopes with 
higher spin constrain more directions in parameter space. For example, out of the seven Xe 
isotopes, only 129 Xe and 131 Xe have non-zero spin and can lead to constraints in the {Og, Og} 
sector. The former of these, 129 Xe has spin 1/2, and consequently its form factors for SD1, 
SD2, and LD receive contributions only from their J = 1 partial waves. E.g. for SD2 (£') 

^•"'V) = ^Eoiisrii^^oiisrib)^)- (34) 

/J + 1 J=0 
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The J = partial wave vanishes for SD2, and since j = 1/2, the sum truncates at J = 1. 
Consequently, viewed as a matrix Fg^' N ■* = C{j x )j^F^' N ^ manifestly has rank 1 and is thus 
sensitive only to one linear combination of c 9 ™' ) and Cg . 131 Xe, on the other hand, has spin 
3/2 and its F^' N ^ is sensitive to two different directions. 

Our evaluations of the nuclear responses of common materials utilized as nuclear targets, 
F, Na, Ge, I and Xe, enable the calculations of expected event rates in this extended param- 
eter space of non-relativistic operators. Here we describe in more details a scan of the four 
dimensional space parametrized by operator coefficients in the Lagrangian, for sectors 
{Oi, O3} and {Os,Og}, as well as for a single operator in each sector but allowing for isospin 
variations. Given that the targets vary in sensitivity along different directions in the pa- 
rameter space, we will show the constraints imposed by various combinations of experiments 
XENON, CDMS, SIMPLE and COUPP. That is, we will plot the intersection of the allowed 
regions of these experiments. For the first two experiments we consider separately high-energy 
and low-energy analyses [10j [HI [12], [13] . 

3.3.1 SI and LSD sector 

2 

Operators which are even under all symmetries, ^— 0\ and O3 can interfere. First, in Fig. 
|2] we show the constraints on this sector (i.e. on c 2 = Y^N=pn c iN + c In)- The various 
curves indicate which experiments are included to obtain a particular constraint. Because 
of its large exposure and sensitivity to SI and LSD responses, the constraint is being driven 
mostly by XENON100, and the combined limit does not change significantly when we subtract 
SIMPLE or COUPP. At low m x < 10 GeV, however, SIMPLE becomes constraining due to 
the kinematics in DM-F nucleus collision. It pushes the constraint from XENON10 and 
CDMS (low threshold) significantly. In order to gauge the importance of interference one can 
compare with the O3 constraint in Fig. [TJ where the constraint on c 2 is more severe, going 
down to 1/ (lOOGeV) 2 for certain DM masses. Thus, interference effects can lessen sensitivity 
by an order of magnitude or more in the cross-section. 

3.3.2 SD2 and LD sector 

Here there are two different sectors: operators q 2 0^ m^O^, and Oq are parity-even and 
comprise one sector, whereas m^Os and Og are a similar, parity-odd sector. In each case, 
the operators containing LD (which is always accompanied by the standard SI coupling) 
are more heavily constrained than those containing SD interactions. Thus, for example, Og 
is less constrained than m^Os, and thus interference between them does not significantly 
loosen the bound on c 2 (except for a small region near DM masses of lOGeV). As is the case 
with the standard SD case, experiments containing elements with unpaired neutrons (such as 
XENON100 and CDMS) need to be combined with those containing unpaired protons (such 
as SIMPLE and COUPP), in order to guarantee greater sensitivity to DM. Indeed, we see an 
almost four orders of magnitude improvement in the constraint when SIMPLE or COUPP 
are combined with XENON or CDMS (Fig. g. 
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Figure 2: Constraints on combinations of experiments for the two interference sectors. 
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Figure 3: A plot showing that for the {Og, Og} sector, XENON (green) , CDMS (brown) and 
SIMPLE (gray), COUPP (purple) are sensitive in orthogonal directions. 



Nevertheless, a comparison between single operator constraint (Fig. [T]) and combined 
limits (Fig. [2]) shows that the combined limit comes mostly from the SD interaction. The 
F M piece in 0$, albeit suppressed by DM velocity v 2 , is still strong enough to impede 8 -0 9 
cancellation. In this sector, interference has a minor effect. 



3.4 Flat directions and DAMA 

Several direct-detection nuclear recoil experiments have observed positive signals that beyond 
their expected background source. While these signals could potentially disappear with better 
background modeling, it is nevertheless interesting to try and interpret the anomalies as 
signals of WIMP scattering. To date, these experiments include DAMA [H] , CRESST [TB] , 
and CoGeNT [16]. Traditional analysis of spin-independent and spin-dependent dark matter 
scattering yield a signal region which is clearly incompatible with strong constraints coming 
from null experiments. Using the effective field theory approach to dark matter analysis, one 
can search through a larger parameter space to see if there is a particular combination of 
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Figure 4: The first two graphs are the low-energy experiment and high-energy experiment 
constraints on {Oi,Os} at m x = 9 GeV. The third graph is the {0$,Og} sector at m x = 10 
GeV, along which direction the DAMA constraint region (red) is almost compatible with null 
experiments. Colors for the other experiments are as in Fig. [3j 
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operators which may resolve the tension between DAMA and the null experiments. We have 
searched this parameter space for regions which are compatible with the null experiments and 
the annual modulation witnessed in the DAMA signal. 

The DAMA detector uses Nal crystals, and the complementary kinematics of the light 
and heavy atoms leads to multiple interpretations of the signal from the detector, resulting 
either from sodium or iodine scattering. The traditional spin-independent 0\ interaction 
thus yields two mass ranges compatible with the DAMA modulation signal. By including 
the full {Oi,03} parameter space, which includes nontrivial interference between operators, 
these mass ranges are expanded to 8-13 GeV and 30-110 GeV. Similarly, for the {0$, O9} 
parameter space, the regions are from 8-17 GeV and from 25 GeV to and past the limit of our 
analysis, 500 GeV. These regions are only useful if they increase the compatibility of DAMA 
with null experiments. 

For the low mass ranges, the preferred region for DAMA is due to sodium scattering. 
Sodium has nuclear spin 3/2. Because each sector has a definite parity, only two of the orbitals 
can contribute to scattering within a given operator sector (Either {Oi,0 3 } or {0 8 ,Og}). 
Thus there are two "flat directions" and two active directions in each sector. This allows 
for a large region of compatibility with DAMA, and by exploiting this freedom it is typically 
possible to fit DAMA with a single given null experiment. 

In the {Oi, O3} parameter space, the strong response of sodium under the $" form factor 
for low ( < 10 GeV) dark matter helps to reconcile DAMA with the null experiments. In the 
low-energy analysis, there is a 1 GeV window centered on 9 GeV which is compatible with 
CDMS, SIMPLE, COUPP, and XENON-10. ( See Fig. [§) However, it is not possible to find 
a region which is compatible with the very strong constraints coming from XENON-100. 

In the {08, Og} parameter space, it is again not possible to find a region that is universally 
compatible. However, the strength of the constraints in this parameter space is heavily 
dependent on the complementary sensitivities of XENON and SIMPLE. If one does not include 
either XENON or SIMPLE, it is possible to find regions which are compatible with the three 
remaining null experiments. This is illustrated in Fig. |4| 

4 Discussion 

Although the results of direct detection searches are usually interpreted only in the standard 
spin-independent (SI) and spin-dependent (SD) scenarios, the full set of possible WIMP- 
nucleon interactions is much richer. Recently, in [TJ, extending previous work by [2J, we 
parameterized the full range of such interactions by constructing the non-relativistic effective 
field theory for dark matter direct detection. In this paper, we have considered analyses 
of the experimental constraints on this parameter space for elastically scattering WIMPs. 
The major qualitative point of these interactions is that besides the SI and SD nuclear re- 
sponses, there are angular-momentum (LD) and angular-momentum-and-spin (LSD) nuclear 
responses, which can differ significantly in their relative strengths at various nuclei from the SI 
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and SD responses. Furthermore, the SD response is a particular combination of two different 
possible spin-dependent responses, one (SD1) that couples to the longitudinal component of 
nucleon spin, and a second (SD2) that couples to the transverse component, and these have 
form factors that differ in detail. 

In section 2 of this paper, we have given a brief review of the effective theory itself and 
how event rates can be calculated for an arbitrary theory using the form factors in [T] for 
the additional nuclear responses that appear. It is our hope that this can serve as a useful 
guide to those wishing to perform an analysis themselves of the constraints on the effective 
theory or some subset of it. In the remainder of the paper, we have analyzed the experimental 
limits on WIMP-nucleon cross-sections that obtain when one considers general points in the 
effective theory parameter space. Our first set of results are the constraints on individual 
operators, shown in Fig. [T], where we have chosen one representative operator for each of the 
possible responses. We have also taken into account the importance of possible interference 
between different responses, which allows cancellations between different responses to occur 
that suppress the dark matter signal at particular isotopes. Our results are presented as the 
constraints on two different subsectors of the effective theory. The first subsector contains 
the SI, SD2, and LSD responses, while the second contains the SI, SD2, and LD responses; 
together, these contain all non- vanishing interference terms between different responses. These 
results are shown in Fig. [2} One relevant question this enables us to address is to what extent 
the limits from combinations of experiments improve upon the limits of individual experiments 
themselves, or in other words, whether experiments can fill each other's "gaps" in sensitivity. 
We find that in the first subsector, the constraints from XENON100 are sufficiently strong that 
not much improvement is gained by combining them with additional experiments. However, 
in the second subsector, we find that, mostly due to the smallness of the proton SD response 
in xenon, significant improvement is made by combining XENON100 with COUPP. 

Finally, we have explored whether the anomaly at DAMA can be interpreted as dark 
matter signal consistent with all other experiments at any point in parameter space. We find 
that, within our analysis assumptions, we rule out any possible elastically scattering solution. 
However, certain highly tuned points in parameter space are ruled out only by a factor of 2 
in cross-section, suggesting that with different analysis assumptions, in particular about the 
halo velocity distribution, a consistent interpretation might be allowed. 
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A Experiments and analysis 



In this section we will briefly describe each direct detection experiment included in our analy- 
sis. To perform a conservative analysis we tag all reported events as WIMP scatterings, except 
where background analysis is specifically mentioned. To calculate the allowed or excluded re- 
gions, either x 2 statistics or maximum gap methods [T7] are used based on the characteristics 
of the data in each experiment. 

A.l DAM A 

In our calculations, the DAMA signal is taken from [TJ]. For sodium scattering events in Nal, 
we take a quenching factor of Q]^ a = 0.3, while for iodine, Qi = 0.085. No channeling effect 
is considered in this calculation. The total degree of freedom of the dataset we use is 17. The 
90% CL. exclusion is therefore the contour with a \ 2 = 24.7. 

A.2 CDMS II 

The CDMS II analysis used germanium detectors searching for WIMP interactions. In the 
data taken between July 2007 and September 2008, with an exposure of 612 kg-days, two 
candidates events at recoil energies of 12.3 keV and 15.5 keV were observed [18]. To calculate 
the 90% CL. exclusion, we use the pmax, or maximum gap method [17], namely setting the 
upper limit of the probability of seeing no events between the energies of the observed events 
and thresholds. 

A. 3 CDMS low energy threshold analysis 

The CDMS collaboration reanalyzed their data taken between October 2006 and September 
2008 to achieve a lower recoil energy threshold [TU] . The search took place within the energy 
region between 2 keV and 10 keV. Events in this region were previously rejected because of 
the high electron recoil background. In their analysis, the nuclear recoils were discriminated 
from electron induced events based on the fact that electrons cause more ionization than 
nuclear recoils of the same energy. Events were identified as coming from a nuclear recoil 
based on their ionization energy and recoil energy, with cuts being determined by calibration 
with 252 Cf sources. 

We take the data as well as the nuclear-recoil acceptance efficiencies provided in the same 
analysis, and combine the events in the eight detectors together. These 400+ events are then 
divided exponentially with respect to the recoil energy into 20 bins. This method of binning is 
chosen as the measured event rate falls quickly with energy. A x 2 statistics is employed with a 
90% constraint at \ 2 = 28.4. In order to account for the unknown background conservatively, 
only bins where the predicted event rate was greater than the observed rate were counted. 
Thus, in bins where the predicted rate is lower than the observation, the difference is posited 
to be due to background. 
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A. 4 XENON10 low energy threshold analysis 

The XENON10 collaboration also performed a low energy study of their XENON10 data, 
down to the nuclear recoil energy of about 1 keVnr (nuclear recoil energy) [12J. The exper- 
iment utilized liquid xenon as the target, which, containing an unpaired neutron, is more 
sensitive to neutron induced scatterings. The detectors recorded scintillation signal (SI) and 
ionization signal (S2). The ratio of the two signals (S2/S1) acted as the criteria of back- 
ground rejection, since electron and nuclear events have different values. However, the energy 
threshold of the scintillation efficiency L e ff for XENON 10 was about 5 keV. To achieve lower 
threshold, they required only observations of the ionization signal and calibrated the events 
using the ionization yield Q y , the ratio of S2 and recoil energy. 

We extract the data from Figure 2 of their analysis, and the ionization yield from Figure 
1. Because of systematic uncertainties, we approximate conservatively, as suggested, Q y = 
4 electrons/keVr below E nr around 30 keVr, such that the newly adopted Q y is a continuous 
function. Beyond the analysis threshold there are five bins of data, containing 41 events, 
among which, the four events with lowest 5*2 values are reported exactly. Since the S2 values 
for the events in the latter bins are unknown, we presume them distributed evenly in each 
bin so as to employ the \ 2 method. 

The resolution of the experiment is taken into account as described in [19J, which presumes 
that the S2 signal follows a Poisson distribution. It is then convolved with the expected event 
rate from 1.4 keV, the cutoff in Q y , to 60 keV, where the event rate becomes negligible. 

A.5 XENON100 

The XENON 100 data acquired for 13 months during 2011 and 2012, with 224.6 days x 34 
kg exposure, exert by far the strongest constraint on the SI scattering cross section in the 
WIMP parameter space [13]. Two events at 7.1 keVnr and 7.8 keVnr were detected, with 
background expectation of (1.0 ± 0.2) events. 

As for the scintillation efficiency L e ff , the parameter relating the nuclear recoil to electron 
equivalent energy and thus the photoelectrons, we directly use the data from the interpolation 
on measurements cited in this study. For nuclear recoil energies below 3 keVnr, there are 
no measurements of L e ff, and so we adopt the logarithmical extrapolation down to zero 
at 1 keVnr. Since the threshold energy in this experiment is 6.6 keVnr, the analysis will be 
relatively insensitive to the particular extrapolation used even after poisson fluctuations of the 
electron signal are taken into account. Also, the resolution is (Jpmt = 0.5 photoelectrons [20] . 
or translated to recoil energy, 1.4 keV. Thus the extrapolation of the scintillation efficiency 
under 3 keV has limited effect on the expected signal events. The uncertainty in L e ff is 
neglected in our analysis. 

The formulae of computing the expected events after considerations of statistics and res- 
olutions is given in [20]. We again impose the 90% C.L. exclusion of our expected event rate 
on the XENON100 results with the maximum gap method. 
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A.6 SIMPLE 2011 



The SIMPLE experiment applies superheated liquid C2CIF5 droplets in bubble chamber-like 
detectors, which is particularly sensitive to spin-dependent scattering due to the fluorine 
content. The first stage of the phase II SIMPLE dark matter search reported a 14.1 kg-d 
measurement in 2010 [21], where 14 events were observed. In 2011 the final results of the 
phase II experiment were presented [22], which re-evaluated the data from the first stage, 
identified 5 out of the previous 14 events as background, and reduced the expected rate to 
0.289 evts/kgd at 90% C.L.. This imposed an even stronger constraint on the cross sections 
than the newer result. A refined efficiency is provided in the same study, which is included 
in our calculations. We integrate over the expected event rate starting from the threshold 
energy of 8 keVnr to reach this value as the 90% exclusion in the parameter space. 

A. 7 COUPP 2012 results 

The COUPP experiment is also sensitive to spin-dependent scattering because of its CF 3 I 
targets. The experiment reported results of the running from September 2010 to August 2011 
with effective exposure of 437.4kg-days [4j. Twenty single nuclear recoil event candidates 
were recorded. In our analysis, the efficiencies of the nuclear targets are taken from |23j . The 
observed events are presumed to follow the Poisson statistics, which are compared with the 
expected event numbers to set a 90% exclusion limit. 

B Electromagnetic moments in the effective theory 

Dark matter interacting with the atomic nuclei through electromagnetic moments is a com- 
pelling and well-motivated special case of the general effective theory, and falls outside of the 
standard spin-independent and spin-dependent interactions. Such electromagnetic moment 
interactions have been considered in the past (the earliest, to our knowledge, being [2"4"]). but 
without the availability of the detailed nuclear response functions. In this appendix, we will 
connect such interactions to our effective theory, which in turn allows one to calculate the 
corresponding momentum- dependent form factors. 

Restricted to the proton and neutron, the EM current is 

Jem = ^ fa + ifoTfc + *f > ( 35 ) 

where p(n) subscripts denote that the operator acts on protons(neutrons). Here, g p = 5.59 
and g n = —3.83. Taking the non-relativistic limit, the leading contributions to are 

Jem = 2mjvl p , (36) 
f EM = K l l p -ig p e^ k q^S k p -ig n e^ k q=S k . (37) 

The first of these is (at this order) boost-invariant (this is just the statement that j EM shifts 
by a term that is suppressed by powers of momenta), whereas the second is not. Under boosts 
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shifting the velocity as v l — > v 1 + /3\ the current j EM transforms according to 

Jem ~^ Jem = Jem ~ ^Jem- (38) 
To make a boost-invariant quantity, we can use P l or K\ which transform according to 

pi _^ p i -2m x p i , K { ->■ JC -2m N /3 i . (39) 
From this we can construct the following boost-invariant: 

— * 

Jem ss Jem - J EM ^ = (-2m^ - S p % otons + ("^ V^)^ . (40) 

— * 

An obvious second choice besides Jem for a boost-invariant would be to use K % instead of P % 
in the above definition. However, this does not produce a new independent invariant, because 
it is related to invariants we have already constructed. Specifically, 

-> K - 

Jem - Jem « = J em + Jem^ • ( 4 1 ) 

ZTYljy 

— » — > 

A similar discussion holds for the electric and magnetic fields E and B themselves. Since q u 
vanishes (at the order we consider), we in fact have 

E = i ~ 2 Jem ■> B = i-^ q x J EM , (42) 
q z q z 

— * — * 

so E is boost-invariant but B is not: 

E ->■ = £, B^B' = B- /3xE. (43) 

— * — * 

This can be remedied as before by constructing a combination of B and E that is invariant: 



B = B-^—ExP = t^qxJ EM . (44) 



Now, it is straightforward to work out the form factor for a general electromagnetic moment 
interaction. Consider the following electromagnetic moments: 

(- S - S -> S 1 A 

UdmB • T^TT + d DME ■ + a DM J EM ■ t^t + -^er 2 D iq- E J , (45) 

which are (boost-invariant generalizations of) a magnetic dipole moment /^dm, an electric 
dipole moment g?dm, an anapole moment Odmj and a charge radius r D , respectively. Applying 
the above results, we see that the leading EFT coefficients for this Lagrangian are simply 

1. Magnetic dipole moment: 



C = 2m xf , DM (^(0 6 -q 2 4 )-^0 5 ) 

W\s x \ t\s x \ /, 

+2m xf i DM (^i{0 % - q 2 A )\ . (46) 

\1 Pxl /neutrons 



protons 
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2. Electric dipole moment: 



2m 



C = 2m x d DM I "','^, 01! 



X I / protons 



(47) 



3. Anapole moment: 



l^xl l^xl / protons 

9n 



+2m x a D M ( 77^7^9 



X I / neutrons 



(48) 



4. Charge radius: 



£ = 2m x er 2 D ( —-m N 0\ 



1 

3 1 



protons 



(49) 



C Form Factors for Effective Theory Operators 



The full set of form factors necessary for the effective operators in eqs. (|3])-([6]) can be written 
in terms of a smaller number of basic independent ones as follows: 
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where 

CUx) = Ajxi3 l + l) . (51) 

All interference terms that are not listed above vanish. (9 2 is omitted as it does not appear 
at leading order from any relativistic interaction. 
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